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Introduction 

The magnetic inverse solution of: 

b = Gq, 

can be transformed into a minimization of: 

[b-Cqf 

with a least squares solution: 
q = [G T G]' 1 G T b 
b = vector of M data channels 
q = vector of N source parameters 
G = matrix of parameter magnetic fields 
with ( M > N) 

Alternatively, the equation can be solved by 
refining an initial solution estimate qo, by a 
gradient descent method [1]: 

q i+ i = qi + a i [G T G]- 1 G T r i (1) 

r i+ i = r f - G(q i+ i - qi) 

The gradient decent technique is commonly 
used to determine an equivalent current dipole 
solution where the source vector incorporates 
both amplitude and location: 

q = [qxq y qzx q y q z q ] 

However, for both source solutions, the 
components of b corresponding to low 
amplitude eigenvalues in spectrum of the matrix 
[G t G] are amplified. This is undesirable, 
especially if b contains noise. As an alternative 
to regularization techniques, the excessive 
sensitivity to small eigenvalue components of b 
can be avoided by performing a source space 
minimization of: 

[G T b - G t G q ] 2 (2) 

with the recursive solution: 

q i+ i = qi + ai[G T G]G T ri (3) 

ri+i = i*i - G(q i+ i - qi) 

oti = i*i T [GG T ] 2 i*i/i*i T [GG T ]V 

This recursive solution amplifies the importance 
of large eigenvalue components of the magnetic 
field rather than the small ones. A sufficiently 
accurate solution is obtained when the number 
of iterative steps is equal to the rank of the 
matrix, G T G, (Eqn. 1). Thus, it is an efficient 
technique if the matrix, G T G, is small. The rank 


of G t G is two at each cortical location if G is 
based on a spherical head model and rank 3 if 
the forward solution incorporates a realistic 
head model such that radial sources generate 
detectable magnetic fields. In addition, 
Equation 3 is utilized in our 2DII imaging 
technique where the rank of G T 2 diiG 2 dii is four 
[2, 3], 

Methods 

The accuracy of this minimization technique 
was tested by using Equation 2 to calculate [q x 
q y q z ] at 16 locations on a 0.6 cm square grid, 
within a MEG sensor array. In Figure 1, a 
spherical volume conductor forward model of 
the head was used to generate 148 channels of 
magnetic field data corresponding to a current 
dipole oriented successively in the x, y, or z 
direction at one grid location. The z orientation 
is almost radial relative to the chosen sphere 
origin. For these three sets of simulated data, 
Equation 3 was used to solve for the solution 
vector, q, at each grid position. In addition, for 
each of the 16 grid points, the fractional 
magnetic field power mismatch of the source 
solution was calculated: 

ERROR =|b-b gridj3t | 2 /|b| 2 . 

and 

goodness of fit = 1 - ERROR. 



Figure 1: left to right,Magnetic field of x, y, and 
z oriented dipole within a small patch of 16 
source imaging locations, (red patch). Black 
dots are at 148 MEG detector locations, (top 
view). 











Finally, imaging software incorporating this 
technique was applied to patient data to 
generate equivalent current dipole (ECD) and 
2DII [1] images of visual evoked cortical 
activity. 

Results 



Figure 2: Goodness of fit contour plots for 
source solutions at 16 grid locations, (black x), 
Active current dipole at black dot with current 
dipole orientation in the x (left), y (center), z 
(right) direction. The Magnetic fields of these 
sources are in Figure 1. The goodness of fit of 
the solution was 100 percent at the grid 
location corresponding to the active source. 
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Table l:The orientation of a unit amplitude 
simulated current dipole source is in first 
column. The x, y, and z amplitudes of the 
corresponding source solution at the active grid 
site are in columns 2-4. The x, y, and z 
components of the radial direction vector of the 
sphere model at this site is included to 
demonstrate that the source solutions are purely 
tangential. Thus, only the small tangential 
component of the z oriented source is 
magnetically visible. 


In Figure 2, the goodness of fit for the solution, 
(Equation 3), is 1.0 at the true source location. 
At each of 16 grid locations, the source 
solutions, q, (Equation 3), were tangential to the 
spherical volume conductor model used to 
generate the simulated data and forward 
solutions. These results are tabulated in Table 1 
for the grid location corresponding to the 
simulated source. The goodness of fit for 16 
source solutions within the 0.6 cm square grid 
ranged from 98.2 to 100 percent at the actual 
source location, Figure 2. For these spherical 
head model results, only 2 iterative steps were 
required to calculate each ECD solution, q. 
Three iterative steps would be required if the 
forward solution incorporated a realistic model 
of the patient's head geometry. 

Figure 3 depicts the results of applying 
Equation 3 to visual evoked field data of a 
patient with a right occipital tumor. The latency 
was 167 milliseconds after checkerboard pattern 
reversal within the left visual field. The three 
dimensional view is derived from the patient 
MRI data. Using the MRI data, the cortical 
gray mater was modeled with a three 
dimensional grid of 2955 grid points and 
Equation 3 applied at each point to determine q 
and the corresponding goodness of fit. At the 
center of the red false color area, the goodness 
of fit was 98.1 percent and the localization error 
estimated to be less than two millimeters. 



Front 


Figure 3: Equivalent current dipole goodness of 
fit metric generated with Equation 3 for left 
visual field data. The adjacent outline of the 
Occipital tumor is visible. 























The 2DII whole brain source activation image 
for this patient data is an average of 32 2DII 
source amplitudes for the time window 167 to 
228 milliseconds after pattern reversal, Figure 
4. The 2DII technique converts the 
underdetermined continuum model of the cortex 
to an iterative solution of a sequence of 
overdetermined source solutions where the rank 
of G T 2 dii G 2 DI 1 is four. 



Figure 4: The average of 32 2DII source 
solutions for 167 to 228 milliseconds after 
checkerboard pattern reversal within the left 
visual field. 

Discussion 

The iterative solution of Equation 2 is a 
conjugate direction minimization technique [4]; 
the sequence of residual magnetic field vector 
changes are orthogonal, ArjArj = 0. Although 
the algorithm operates on the residual field, r i? 
the final solution can be represented as: 

q = E[3i[G T G] 2i ' 2 G T b, i >2. 

This form of the solution can also be obtained 
by using the characteristic equation of the 
matrix [G T G] to construct a series expansion for 
[G t G] _ 1 in powers of [G T G]. 


For these single current dipole results, the set of 
equations required for the series expansion is: 

[G t G] _1 (P 3 [G t G] 3 + |3 2 [G T G] 2 +Pi[G t G] +|3o} = 0 
... { ... + ... + ... + ... } =0 
[G T G] 6 {P 3 [G T G] 3 + p^Gf+p^G] +p 0 } = 0 

p k are characteristic equation coefficients. 

However, the iterative construction of the 
solution using Equation 3 is efficient and easy 
to adapt to matrices of different rank. In 
addition when applied to real MEG data the 
number of terms in the solution can be truncated 
when the residual field is primarily noise. The 
application of this technique to calculate 2DII 
source solutions and equivalent current dipole 
source amplitudes demonstrates the efficiency 
of this technique when the matrix, [G T G], is 
small 
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